Variational principle for the Wheeler-Feynman electrodynamics 



Jayme De Luca* 
Universidade Federal de Sao Carlos, 

Departamento de Fisica 
Rodovia Washington Luis, km 235 
Caixa Postal 676, Sao Carlos, Sao Paulo 13565-905 
Brazil 

(Dated: January 8, 2009) 

We adapt the formally-defined Fokker action into a variational principle for the electromagnetic 
two-body problem. We introduce properly defined boundary conditions to construct a Poincare- 
invariant-action-functional of a finite orbital segment into the reals. The boundary conditions for the 
variational principle are an endpoint along each trajectory plus the respective segment of trajectory 
for the other particle inside the lightcone of each endpoint. We show that the conditions for an 
extremum of our functional are the mixed-type-neutral-equations with implicit state-dependent- 
delay of the electromagnetic-two-body problem. We put the functional on a natural Banach space 
and show that the functional is Frechet-differentiable. We develop a method to calculate the second 
variation for orbital perturbations in general and in particular about circular orbits of large 
enough radii. We prove that our functional has a local minimum at circular orbits of large enough 
radii, at variance with the limiting Kepler action that has a minimum at circular orbits of arbitrary 
radii. Our results suggest a bifurcation at some 0(1) radius below which the circular orbits become 
saddle-point extrema. We give a precise definition for the distributional-like integrals of the Fokker 
action and discuss a generalization to a Sobolev space of trajectories where the equations of 
motion are satisfied almost everywhere. Last, we discuss the existence of solutions for the state- 
dependent delay equations with slightly perturbated arcs of circle as the boundary conditions and 
the possibility of nontrivial solenoidal orbits. 



INTRODUCTION 



We construct a variational principle for the electromag- 
netic two-body problem with finite integration limits. 
Unlike the Fokker action that involves an infinite integra- 
tion and has a formal meaning only[l~3], our Poincare- 
invariant functional maps a finite segment of trajectory 
into a finite real number. Our variational principle uses 
boundary conditions consisting of (i) the initial point Oa 
for the trajectory of particle 1 plus the segment of tra- 
jectory of particle 2 inside the lightcone of Oa, and (ii) 
the endpoint Lb for the trajectory of particle 2 plus the 
segment of trajectory of particle 1 inside the lightcone 
of Lb. For trajectories respecting the above bound- 
aries we show that the conditions for an extremum of 
our functional are the two-body equations of motion of 
the Wheeler-Feynman electrodynamics [3] . Our first func- 
tional is the natural generalization of the Fokker action 
and it can not be defined for trajectories travelling faster 
than light (superluminal). We construct a norm for the 
linear space of orbits satisfying the above bound- 
aries and show that our functional is Frechet differen- 
tiable at subluminal orbits. In order to obtain a func- 
tional defined everywhere on a natural Banach space of 
orbits we give up the parametrization-independence 
and construct a second generalized functional. The elec- 
tromagnetic equations of motion follow from the con- 
ditions for an extremum of our functionals in the sub- 
space of orbital variations. The extremum condi- 
tions are parametrization-independent for the Fokker- like 



functional, while for the generalized functional the con- 
ditions require the parameter to be proper-time because 
of a conservation law that separates the extremal orbits 
in three classes. The domain of our second functional is a 
Banach space, and moreover along its extremal orbits the 
particle trajectories turn out to be of three possible types 
(a) subluminal trajectories travelling slower than light ev- 
erywhere, (b) luminal trajectories travelling at the speed 
of light everywhere or (c) superluminal trajectories trav- 
elling faster than light everywhere. We calculate the sec- 
ond variation of the action about general orbits for or- 
bital variations satisfying the above boundary conditions 
and in particular about the Schoenberg-Schild-circular- 
orbit-solutions of a large radius [9, 10]. We prove that the 
second variation is positive-definite about circular orbits 
of large enough radii, so that circular orbits are local min- 
ima of our functionals. Our results suggest a bifurcation 
at some 0(1) radius below which circular orbits become 
saddle-point extrema, at variance with the Kepler action 
for which circular orbits of arbitrary radii are minima[12]. 
We discuss a use of the variational principle to solve the 
neutral-delay equations of the electromagnetic two-body 
problem as a boundary-value problem with a variational 
integrator [4]. We discuss a generalization to a Sobolev 
space -ffg of trajectories where the equations of motion 
are valid almost everywhere and the existence of solu- 
tions with slightly perturbed circular boundaries. Last, 
we discuss the physics of the Fokker action and the exis- 
tence of nontrivial solenoidal orbits. 

The Fokker action functional is a synthetic principle of 
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electrodynamics discovered in the early 20th century [1, 2] 
and used in 1945 by Wheeler and Feynman[3] to con- 
struct an electrodynamics of point charges. The Wheeler- 
Fcynman electrodynamics is an alternative description 
of classical electromagnetism that avoids the notion of 
field to describe the classical laws of Gauss, Faraday, 
Ampere, and Biot-Savart [3, 5]. The theory describes 
point charges interacting in pairs via the half-retarded 
plus half-advanced solutions of Maxwell's equations for 
the fields [6]. Here wc avoid the popular name action- 
al- a- distance electrodynamics because it can suggest 
action-at-the-same-time connecting spatially-separated 
points, while the Wheeler- Feynman theory involves only 
Einstein-local interactions along lightcones. Among the 
existing versions of electrodynamics of point charges [7], 
the selling points of the Wheeler-Feynman theory are 
(i) The point-charge- limit is regular, i.e., a spherical 
charge distribution of a small radius does not make a 
force on itself and its mass is not renormalized and (ii) 
The theory reduces to the usual Dirac electrodynamics 
with retarded-only interactions [8] when the far fields van- 
ish asymptotically, a condition named the absorber hy- 
pothesis in Refs. [3]. The equations for two-body mo- 
tion of the Wheeler-Feynman theory are state-dependent 
neutral-delay equations and little is known about their 
solutions, besides the existence of a one-parameter fam- 
ily of circular-orbit solutions[9, 10]. An existence result 
was proved in Ref. [17] for the two-body problem with 
equal charges (repulsive interaction) and initial condition 
restricted to colinear orbits of large separations, a case 
where the equations of motion are no longer neutral but 
rather delay-only. References [18, 19] considered satisfy- 
ing the state-dependent neutral-delay equations almost- 
everywhcre. In Ref. [20] the equations of motion were 
expressed as an algebraic-differential system by solving 
for the most advanced accelerations, an approach also 
used in Ref. [21] to prove wcllposcdncss and existence 
for C°° initial data consisting of maximal independent 
past segments. The initial conditions consisting of maxi- 
mal independent history segments developed in Ref. [21] 
are different from the initial conditions used in our vari- 
ational method, which combine future and past data. 
Last, the simpler delay-only state-dependent two-body 
equations with initial condition restricted to colinear or- 
bits were studied numerically in Ref. [25] for the case of 
repulsive interaction and in Ref. [26] for the case of at- 
tractive interaction (opposite charges). This paper aims 
to introduce the problem for a mathematical audience. 
In the introduction we start from the naive and formal 
language of physics, posing the problem first at an intu- 
itive level. In the later sections we make an attempt to 
proceed with rigor and precise definitions by presenting 
the results in the form of theorems. 

The paper is divided as follows: In section 1 we give a 

crash review of Minkowski spaces and put in one place the 
ingredients later used to show that the equations of mo- 



tion separate the orbits in three invariant classes and to 
construct an action defined on a Banach space. In section 
2 we introduce the finite action and the boundary condi- 
tions. Wc construct a norm enforcing the property that 
perturbations with a small norm of subluminal orbits 
yield subluminal orbits. For such norm the functional 
is Frcchct-difFcrcntiablc along subluminal orbits. In or- 
der to obtain a functional defined on a complete normed 
space we relax the parametrization independence require- 
ment and construct a second functional that can be ex- 
tended to all types of orbits of a natural Banach space. In 
this section we discuss the advantages of using the varia- 
tional method as an alternative to the state-dependent 
neutral-delay equations of motion as far as numerical 
stability. In section 3 we give a method to calculate 
the second variation about arbitrary solutions for or- 
bital perturbations. In particular we calculate the second 
variation about low-velocity-circular-orbit-solutions. We 
show that the quadratic form of second variation about 
low- velocity-circular-orbits is positive-definite if the cir- 
cular radius is large enough. In this section we develop 
the idea of a sewing grid which appears naturally in the 
integration of the quadratic form with delay and is use- 
ful for the numerical analysis of state-dependent delay 
problems. In section 4 we put the discussions and the 
conclusion. We discuss the variational method as a tool 
to investigate solenoidal and other types of orbits of the 
electromagnetic two-body problem. We also discuss the 
variational problem with slightly perturbed circular-orbit 
boundary data. Last, in the appendix we review the 
physics of the Fokker action and the conserved momenta 
of Noethcr's theorem. Wc discuss the nontrivial possibil- 
ity of solenoidal orbits with both particles gyrating near 
the speed of light with finite and small momenta. 

PRELIMINARIES AND DEFINITIONS 

We start by explaining the natural coordinatization for 

Lorcntz-invariant dynamics, i.e., the Lorcntz four-space 
L'^ attached to an inertial frame by Einstein synchro- 
nization of clocks (the L in L'^ stands for Lorentz). A 
point in L'^ is defined by a time t and a spatial position 
r in the inertial frame, = {t,f), henceforth called 
the time component t and the three-vector spatial com- 
ponent r. The index /i belongs to (1, 2, 3, 4), with /i = 1 
denoting the time-component and /x = 2,3,4 denoting 
the spatial components. From any Minkowski vector 
a'^ = (ao,a) wc define its dual vector bya^ = (ao,— a). 
The Minkowski scalar product is a bilinear product de- 
fined as the usual scalar product on between the first 
vector and the second vector's dual (or vice- versa), i.e., 
(a • 6) = ai^i — 0262 — a3&3 — aibi . This definition gives 
only a pseudo-scalar bilinear product, and the pseudo- 
norm |a|^ = (a ■ a) induced by the Minkowski product is 
sensible only for time-like vectors, i.e., when (a • a) > 0. 
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The Minkowski product divides the vectors of L"* in three 
classes : (i) if (a • a) > the vector is called time- 
like (for example the four-velocity along a subluminal 
orbit), (ii) if (a • a) < the vector is called space-like 
(for example the four-acceleration of a subluminal orbit) 
and last (iii) if (a • a) = the vector is called a null- 
vector or light-like. The four-vectors a and b arc said 
to be orthogonal if {a ■ b) = 0. The properties of the 
Lorentz group and the Minkowski product are discussed 
in Ref. [23], of which wc list a few:-(a) Two orthog- 
onal light-like vectors are necessarily multiples of each 
other because (a • a) = {b ■ b) = {a ■ b) = implies 
the Cartesian product of the three-vector components 
satisfies fa'^a = ||?a||||n)|| (double bars denoting the 
Euclidean modulus of the three- vector), (b) All vectors 
orthogonal to a time-like vector are space-like and form a 
three-dimensional space-like subspace. (c) Given a time- 
like four- vector z and an arbitrary four-vector x there is 
a unique decomposition x = y + a.z, where y is space- 
like and a a real scalar, and (d) any orthogonal basis 
fori^ must contain one time-like four-vector and 3 space- 
like four- vectors [23]. (e) For time- like vectors the invari- 
ant reverse-Schwartz-inequality holds for the Minkowski 
product, i.e., (a • b)^ > (a • a)(6 • 6), equality holding 
iff the vectors are parallel, and last (f) For a time-like 
and a space-like vectors, the reverse Schwartz holds with- 
out the equal-sign case, i.e., (a • 6)^ > (a • a)(b • b) [23]. 
The positivity of the Minkowski product (a ■ a) for the 
four-velocity in arbitrary parametrization is the physi- 
cal condition that the particle travels slower than light. 
The four-velocity is light-like in the limit situation when 
the particle travels at the speed of light. The Minkowski 
scalar product is left invariant by Lorentz transforma- 
tions and it is useful to express the equations of motion 
and the action functional in a form explicitly equivari- 
ant under the Lorentz group. The last three components 
of a Minkowski vector form a spatial three- vector usually 
treated differently from the first component, and the var- 
ious norms used in this papers are henceforth denoted as 
follows: (I) the Minkowski norm is indicated with single 
bars, i.e., |a|, (II) the Euclidean M'^norm of the spatial 
three-vector part is indicated by double bars, i.e., ||ra||. 
We also use double bars to indicate the absolute value 
of a real number and (III) the Euclidean norm of a 
four-vector is indicated by double bars with sub-index 4, 
i.e., ||a||4 and last (IV) The norm defined on our func- 
tional linear space of trajectories of section 3 is denoted 
by |bi|jv(xi)- 

To abbreviate the notation, we henceforth drop the 4- 
index label and keep only a lower index j G (1, 2) to iden- 
tify each particle of the two-body problem, e.g., j = 1 de- 
notes electronic four-vector quantities and j = 2 denotes 
protonic four-vector quantities. For subluminal orbits it 
is convenient to express the equations of motion in terms 
of a Lorentz-invariant parameter defined by the squared- 
Minkowski-norm of the infinitesimal displacement vector 
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{dnf = {dtif - {dxif - {dvif - {dzif > 0. (1) 

The left-hand side of Eq. (1) is positive for subluminal 
orbits, zero for luminal orbits and negative for superlu- 
minal orbits. The parameter Tj defined by Eq. (1) is 
called the proper-time and it is a property of each parti- 
cle's trajectory, the usual parametrization by arc-length 
of differential geometry. 

Next we introduce the naive Fokker action in the above 
defined Lorentz four-space using a normalized unit 
system where the speed of light is c = 1 and the elec- 
tron and the proton have mass and charge mi = 1 and 
ei = — 1 and m2 = 1824 and 62 = 1 respectively. Let 
the trajectory of each particle in be a differentiable 
function Xi{\i) : M ^ of a parameter \i with i = 1,2 
indicating respectively the electron and the proton tra- 
jectories. The Fokker action[l, 2] is defined in the orig- 
inal literature by a formal integration along the whole 
trajectories as 

S = — J mi^/xi ■ xidXi — J m2^/x2 ■ X2d\2 

+ j j 5{\xi - X2f)xi ■ X2dXidX2, (2) 

where overdot denotes derivative respect to the param- 
eter of each trajectory. Action (2) is formally indepen- 
dent of the parametrizations, a geometric property easily 
checked by changing the parameter of each trajectory 
with the chain rule. The peculiar last integral of the 
right-hand-side of Eq. (2) involves the composition of 
the Dirac delta-function S{x) with the real function 

d{XiX2)=\xi{Xi)-X2{X2)\\ (3) 

where single bars stand for the Minkowski norm of the 
four-separation X12 = {xi — X2)- The peculiar combina- 
tion appearing in Eq. (2) comes from the Green's func- 
tion of Maxwell's equations and becomes ill-defined along 
trajectories or in a Sobolev space. Here we give a pre- 
cise definition for the right-hand-side of Eq. (2), and in 
Section 3 we define a consistent derivative for such opera- 
tion before evaluating the second variation, thus avoiding 
the loose derivatives of the Dirac delta-function. Condi- 
tion (1) is sufficient for the separation d{Xi, A2) of Eq. (3) 
to have precisely two zeros for each fixed Ai along a 
trajectory [21]. In Ref. [21] it is proved that along a sub- 
luminal orbit light captures the slower moving particle 
once in the past and once in the future. The integration 
on the last term of the right-hand-side of Eq. (2) gives 
a nonzero contribution at each zero (Ai, A2) of Eq. (3). 
At a given Ai the condition 

|a;i(Ai)-a;2(A2)|2 =0, (4) 

can be solved for the time-component t2(A2) of vector 
X2{X2), yielding a retarded time and an advanced time, 
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each defined implicitly by 



i2(A2)=ti(Ai)T||r2(A2)-ri(Ai)||. 



(5) 



where double bars stand for the Euclidean norm of the 
spatial separation r2{^2) — n.(Ai). Either one of the 
equivalent Eqs. (4) or (5) are henceforth called the light- 
cone condition. Equation (5) is an implicit condition 
for A2 because A2 appears on both sides as an unknown 
argument. Condition (4) is symmetric on particle tra- 
jectories, so that the lightcone condition for the protonic 
trajectory is still Eq. (5), as obtained by either rear- 
ranging Eq. (5) to isolate ti on the left-hand-side or by 
exchanging the indices 1 and 2 of Eq. (5). In the fol- 
lowing we assume the orbital parameters are restricted 
to the intervals [Lai, L^f] for a = 1,2, as defined in the 
next section. Let the zeros of (i(Ai, A2) for A2 G [L21, L2f] 
and a fixed Ai e [Lii,Lip] be (Ai, A2''^(Ai)). Definition 
1 : We henceforth define the integral involving the Dirac 
delta- function composed with d(Ai, A2) by 



J(d(Ai,A))/(Ai,A)dA 



L21 
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/(Ai,A<^)) 
|^(Ai, A^-'^) 



(6) 



where ||^^(Ai, A2"''')|| is the absolute value of the partial 

derivative of d{Xi,X2) evaluated at each zero (Ai,A2''^) 
of rf(Ai,A2). The summation on the right-hand-side of 
Eq. (6) includes all the zeros of the lightcone condition 
inside [L2i,L2f]- Once the separation d(Ai,A2) is com- 
pletely symmetric on particle quantities, definition (6) 
has a symmetric definition as follows;- Definition 2: We 
henceforth define the integral over Ai G [Lu, Lip] involv- 
ing the Dirac delta- function composed with d(Ai, A2) by 



Lip 

I 



5(d(A,A2))/(A,A2)rfA = E 



where 1 1 £^ {X^^^ > A2) 1 1 is the absolute value of the partial 

derivative of (i(Ai, A2) evaluated at each zero (A^'^^ A2) of 

(^(Ai, A2) for a fixed A2 G [L2I,L2f] and Ai G [Lu, Lip]. 
For subluminal orbits the interval [Lki,LkF] can include 
at the most the two zeros proved in Ref. [21], while 
for superluminal orbits there can be several zeros inside 
[Lkj, Li;f], or even none. If no zero exists in the integra- 
tion interval the right-hand side of either Eqs. (6) or (7) 
is defined to be zero. Definitions (6) and (7) are moti- 
vated by the evaluation of the respective left-hand sides 
of Eqs. (6) and (7) using the Dirac delta-function with 
a C°° separation d{Xi,X2) and changing variables using 
the absolute value of the Jacobian of the local coordi- 
nate change near each zero. Here we avoid distributional 
operations with the Dirac delta-function and henceforth 
take Eqs. (6) and (7) as defining a functional of (i(Ai, A2), 



/(Af\A2) 



,A2 



(7) 



/(Ai, A2) and the intervals [Lki, Lkp] into the reals. It is 
further useful to define the function A{xj) : L'^ ^ L'^ by 



Ak{xj{X. ))= j 5{[xk{Xk) 



Xj{Xj)[^)xk{Xk)dXk, 



(8) 

where (fc, j) = (1, 2) or (2, 1) and the integration on the 
right-hand side of (8) is defined either by Eq. (6) or Eq. 
(7) . The vector fimction defined by Eq. (8) is often called 
the vector-potential in physics. Assuming integral (8) to 
exist for both {k,j) = (1,2) and (2,1), the interaction 
double-integral of the right-hand-side of Eq. (2) can be 
expressed by 



I = J Ak{xj) ■ x.dXj, 

Lji 



(9) 



with either (fc, 7) = (1, 2) or (2, 1). Using either definition 
(6) or (7) we can express the interaction term (9) in the 
two equivalent forms 



I = 



L21 k 



-^i'^l :A2 



(10) 



I d\ 



I rfAiE 



Lii 



XijXi) ■ X2{X^i 
W 9d {\ 



It is instructive to check the equivalence of formulas (10) 
by changing the integration variable from A2 to Ai about 
each zero. Condition (4) defines Ai as a function of A2 by 
the the implicit function theorem and the Jacobian of the 
coordinate change transforms the first line of Eq.(lO) into 
the second line of Eq.(lO). Last, to express the Jacobian 



in the usual form of physics textbooks we define 



-2dX2^^"^'^ 

[xiiXi)-X2{Xt)]-X2{Xt). 



(11) 



For subluminal orbits X2{Xf) is a time-like vector with 
a positive time-velocity, and once a;i(Ai) — X2(A^) is a 
null- vector, condition 11 defines a positive J^^ on the 
retarded lightcone and a negative J^^ on the advanced 

cone. For superluminal orbits can have any sign in 
either lightcone, so that it is best to keep the moduli in 
the denominators of (6). 



ACTION WITH BOUNDS 

The guiding principle to construct an action functional 
is that the extremum condition should generate the elec- 
tromagnetic equations of motion [1-3]. In the following 
we start from the naive Fokker action (2) and explain 
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how to restrict the integration to suitable finite segments 
of trajectory, using particle-time par ametrization just for 
simplicity of the exposition. The original works [1-3] ex- 
tended the integration of (2) from plus to minus infinity 
as a simple solution to include the needed future or past 
of the other particle at endpoints. There is no reason 
to assume such integration should converge, so that the 
infinite integral (2) has a formal-only meaning[15]. More- 
over, the Fokker action yields the electromagnetic equa- 
tions of motion[l 3] only if the condition of extremum is 
enforced formally with trajectory variations of compact 
support. Here we avoid the shortcomings of a formal- 
only action and give instead a finite-valued functional. 
The boundary conditions can be restricted to a point and 
a segment along each trajectory in a way that the future 
and the past lightcone points exist everywhere along both 
trajectories, as follows; Let the initial point of trajectory 

1 be point Oa at = and the endpoint of trajectory 

2 be point at i2 = T2 as illustrated in Fig. 3.1. The 
trajectory of particle 1 to be varied extends from Oa to 
point L~ at ti = Ti where trajectory 1 intersects the 
past lightcone of Lb { indicated in green in Fig.l). The 
future history of particle 1 is needed from point L~ to 
point Z/+ at ti = > Ti where trajectory 1 intersects 
the advanced lightcone of Lb (the red portion of the up- 
per trajectory of Fig. 3.1). The past history of particle 
2 is needed from point 0~ at t2 = where trajectory 
2 intersects the past lightcone of Oa up to point 0+at 
t2 = AJ < T2 where trajectory 2 intersects the future 
lightcone of Oa (also indicated in red in Fig.l). The tra- 
jectory of particle 2 to be varied goes from 0+ to Lb- 
The combination of the initial point Oa along trajectory 
1 and the final point Lb along trajectory 2, plus the re- 
spective segments of trajectory inside the lightcones of 
these endpoints is henceforth called exchange-of-history 
boundary conditions (EHBCs) as indicated in red in Fig. 
3.1. Our construction is Lorentz-invariant because light- 
cones are Lorentz-invariant objects. The construction is 
unique up to a time-reversed construction using an end- 
point along trajectory 1 and an initial point along tra- 
jectory 2 plus the history segments inside the respective 
lightcones. For orbits the EHBCs complete the trajec- 
tories in such a way that any point along each trajectory 
has the two lightcone roots inside the evaluation interval 
for either one of the interaction formulas (10). 

The restrictions for the EHBCs histories are; (a) it 
must be possible to travel from the initial point to the fi- 
nal point of each trajectory at a speed lesser (or equal at 
the most) than light, and (b) The minim,ally short condi- 
tion that trajectory 1 must intersect the future lightcone 
of 0+ before arriving at endpoint L~ (at time t\ =T{). 
In this way the past history of particle 2 does not inter- 
act with the future history of particle 1. Beyond that 
the variational method can be postulated with otherwise 
arbitrary histories. The advantage of solving the state- 
dependent delay equations using the variational method 



with the EHBCs is the numerical stability:-Once the 
equations of motion are time-reversible the stable and 
unstable manifolds exist in pairs, so that tying both ends 
down with the EHBCs avoids the orbit to diverge either 
in the future direction along the unstable manifold or in 
the past direction along the stable manifold. Since the 
maximum spatial velocity is c = 1, the spatial position 
of particle 1 is bounded by a sphere of radius Ti cen- 
tered at point Oa, while the spatial position of particle 
2 is bounded by a sphere of radius (T2 — Aj) centered 
at 0+. Therefore the subluminal trajectories satisfying 
the EHBCs are spatially bounded and there are no run- 
away orbits satisfying the EHBCs during the optimiza- 
tion. The interaction formula of Eq. (10) needs the posi- 
tion in lightcone along the other trajectory, which is nat- 
urally approximated numerically using the trapezoidal 
rule with an integration grid consisting of the union of 
sewing chains defined as follows; — (i) A forward sewing 
chain is a set of conseciitivc points in lightcone starting 
from an arbitrary point on the boundary segment from 
O- to 0+ (as illustrated in Fig. 4.1). The chain goes 
up to the corresponding point in future lightcone along 
trajectory 1 and back down and up until the last point 
along the boundary segment of trajectory 1 from L~ to 
L+ and (ii) A backward sewing chain is a set of consec- 
utive points in lightcone starting from any point on the 
boundary segment from L~ and 1/+ of the trajectory 1 
(as illustrated in Fig. 4.1). The sewing chain proceeds 
to the corresponding point in past lightcone along tra- 
jectory 2 and back down and up until the last backward 
point on the boundary segment from 0_ to . It is im- 
portant to include in the sewing grid the forward chain 
starting from Oa and the backward chain starting from 
Lb because these chains separate boundary data from 
orbital data. Notice that a sewing chain starts from a 
point along one trajectory and ends with a point along 
the other trajectory, so that each chain defines the same 
number of points along each orbit. 

For arbitrary trajectory variations satisfying the 
EHBCs the linearized functional variation is a sum of 
the linear variations along the two special cases;- (i) one 
fixes trajectory 2 while varying trajectory 1 arbitrarily 
and (ii) one fixes trajectory 1 while varying trajectory 2 
arbitrarily, so that it suffices to study problems (i) and 
(ii). In the following we study (i) using particle-time 
parametrization, for which we integrate over t2 in the 
double integral of action (2) with the help of Eq. (6). 
The half-Jacobian needed for Eq. (6) is a case of Eq. 
(11) with the choice of parameter A2 = t2, i.e., 

Jt, = [Xl{ti)-X2{t2)]-X2{t2). (12) 

The dot over X2 in Eq. (12) denotes derivative respect 

to particle-time. Using (6) to integrate over A2 = ^2 
inside the double-integral on the right-hand side of Eq. 
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(2) yields 



pTi /" J 2 

S = —m\\Jx\ ■ Xidti + / — TO2\/*^2 • X2dt 
Jo Ja+ 

'■^1 ±1 ■ ±2+ 



rT2 



+ 



2II4II 



■dti 



A+ . 

''i Xi ■ X2- 



21 1 J"^ 



■dti, 



(13) 



t2l 



where the superscripts ± on J^^ indicate evaluation on 
the advanced/retarded light-cone of particle 1, respec- 
tively. Notice that action (13) is defined only for sub- 
luminal and luminal orbits. To evaluate the functional 
derivative of S with respect to variations of trajectory 1 
we can drop the last term of the first line on the right- 
hand side of Eq. (13), which is independent of trajectory 
1. Last, the integration over the future history of particle 
1 extending from t\ = Ti to ti = is left invariant by 
a variation of trajectory 1 respecting the EHBCs, so that 
we can replace the upper limit of the last integral on the 
right-hand side of Eq. (13) by Ti, yielding an integra- 
tion over ti G [0, Ti] of a Lagrangian function Li{xi,xi) 
defined as, 



Si = 



j Li{xi,xi)dti 



Ti 

fr r. , Ai-A2+ , a;i • ^2- 

= / {-"mw/xx -Xx-V 



2||4I 



21 U, 



(14) 

\dtx. 



ti I 



In Eq. (14) the advanced/retarded position and velocity 
of particle 2, indicated by ±, are evaluated with the fixed 
trajectory of particle 2 at the advanced/retarded points 
defined by the roots ^2(^1) of Eq. (4), which are implicit 
fimctions of the updated trajectory of particle 1. Notice 
that even though the trajectory of particle 2 is fixed, 
the corresponding lightcone points move away from 
along the fixed trajectory 2 as we vary the trajectory 
1. The gradient of respect to the four-position aji is 
obtained relating the differential dx\ along trajectory 1 
to the differential 
condition (4), i.e.. 



to the differential dt^ via the derivative of the implicit 



-2jfdtf + 2xf2 ■ dxi = 0, 



(15) 



where J^is defined by Eq. (12) and xf2 = xi{ti) — 
X2{tf). Therefore the derivative of tfixi) respect to xi 
along the fixed orbit of particle 2 is 



x 



12 



(16) 



Next we construct a linear space consisting of the 
orbital neighborhood of any subluminal orbit xi sat- 
isfying the EHBCs. Without loss of generality we operate 
with perturbations of trajectory 1 only, defined as type 



(i) in the paragraph above Eq. (12). Definition 3: For a 
subluminal orbit X\ satisfying the EHBCs we define 
the linear space N^^\xi) as the set of all trajectories 
defined by a perturbation function 61 : [0, \if\ — * -^^ , 



I.e., 



■ill 



Xi + bi, 
xi + bi, 



(17) 



where 61 vanishes at the endpoints in accordance with 
the EHBCs, i.e., 



biiXi 
bi{Xi 



0) = 0, 
Aif) = 0. 



(18) 



Notice that the EHBCs forbid orbital perturbations 
61 (Ai) with a monotonically increasing time-component 
because condition (18) is impossible for a monotoni- 
cally increasing time-component. The usual norm for 
the space of functions is given by || = sup||6i||4 -|- 
sup||6i||4 and because of conditions (18) it turns out 
that sup||bi||4 < |Aii?| sup II61II4 for < Ai < Xip as 
can be shown using either one of conditions (18). For 
example using 61 (Ai = Aif) = we have 



ybi{x)dx, 

Ai 



(19) 



so that sup||6i||4 < |Ai^| sup II61II4 for < Ai < Ai^. 
Therefore we can drop the sup||6i||4 term of the norm 
and henceforth our norm is simply defined by the sup of 
the Euclidean norm of bi, i.e., |bi|7v(2;i) ^ sup||bi||4. 
Notice that whenever bi = the condition sup | |6i 1 14 = 
plus the endpoint condition (18) selects the single con- 
stant element bi = 0, so that sup | |bi 1 14 defines a norm on 
the linear space of functions 61 : [0, Xip] satisfy- 
ing the EHBCs. The linear space N'^^\x\) can be shown 
to be a Banach space with this norm in the usual way. 
Proposition 1 : Subluminal orbits (fi(Ai), fi(Ai)) have 
small neighborhoods in N''^\xi) containing only sublu- 
minal orbits. To show it we define the local Cartesian 
velocity respect to particle-time by 



v{Xi)-. 



_ {dn/dXi) 
' {dtx/dXi) ' 



(20) 



a three-vector function of Ai with Euclidean norm lesser 
than one by condition (1). Along a subluminal orbit con- 
dition (1) is positive on the compact set [0, Aif], so that 
(l-||wi|P) > (5 > on [0, Xip\. For any orbit a;i(Ai) sat- 
isfying Eq.(l) we can further define hi = (dti/dXi) > 
and express the velocity ii by 



xi = {hi,hivi). 



(21) 
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Given a perturbation bi = (b, bvi,) and a subluminal orbit 
xi = {hi, hivi), substitution of (17) into (1) yields 

(ii+6i)2 = {hi+bf ^\\hivi+bvb\\^ = 

hiii-\\vi\\^)+b'{i-\\v,r) 

+2hib{l - vi ■ Vb), (22) 

Since the positivity of Eq. (22) is independent of 
monotonic reparametrizations, in the following we use 

paramctrization by the time-component of so that 



hi 



1. The norm |5i 



Nixi) 



dominates the absolute 



value of the time-velocity perturbation b defined above 
Eq. (22), i.e., |6i|iv(a;i) — 1^1' ^° ^^^^ ^^"^^ hmit |&| 
by choosing 5/4 > |bi|Ar(2;j) > |^|. Equation (22) with 
(5/4 > \b\ and hi = 1 shows that the perturbed element 
is subluminal for small enough |6i|jv(xi)) ^^^^ sublu- 
minal orbits have small neighborhoods containing only 
subluminal orbits. 

Next we define the Prechet derivative of action (13) 
about a subluminal orbit Xi: Let 5*1(61,61) : N^^''{xi) 
— » M be defined by substituting (17) into Eq. (14) and 



expanding to linear order for small |6 



l\N{xi)- 



The lin- 



ear expansion of 6*1(61,61) in terms of 61 and 61 in- 



volves integrals controlled by an 0(|6i 



N{xi 



the 



above defined subliminal neighborhood |6i|^(j.j-) < (5/4 
because the Euclidean norm ||6i||4 is also boimded by 
|Aii?| sup II61II4 as explained above Eq. (19). The hnear 
expansion of Si is already the desired Frechet derivative, 



I.e., 



5Si 



J ^dxi 



31 

■bi + —^-bi\dXi, (23) 

OXl 



Even though the functional is already Frechet differen- 
tiable in N'''^\xi), the electromagnetic equations require 
at least a orbit, as follows ; For a orbit Xi the 
second term on the right-hand side of Eq. (23) can be 
further integrated by parts using (18) to yield a term 
linear in 61, so that 5 Si becomes 



5Si 



-I 



Gi ■ bidXi, 



(24) 



with 



dXi dxi dxi ' 



(25) 



where Li is defined by Eq. (14) and Gi E 
is defined only along any orbit of the natural 
neighborhood N^'^\xi). Notice that Eq. (14) is inde- 
pendent of the paramctrization and the expression of Li 
in terms of Ai is obtained simply by replacing ti with Ai 



in Eq. (14). Expressing Eq. (14) as a function of Ai and 
evaluating Gi with Eq. (25) yields 



Gi 



d 

— (mi 
aAi y/xi 



Xi 



■ Xi 



d xi- ±2+ 

dxi^mtw 



_ X2+ 

II A2 1 

Xi ■ ±2- 



X2- 



21 U: 



r) (26) 



2IIJ. 



A2I 



where the dot over Xi denotes derivative respect to A^ for 
i = 1,2 and A2 is the arbitrary parameter of trajectory 
2. The condition for an extremum that follows from Eq. 
(24) is Gi = (plus the symmetric condition G2 = ob- 
tained by varying trajectory 2). Notice that Gi must be 
zero only in the open interval (0, Xip) because the inte- 
grand of Eq. (24) vanishes at the boundaries with 61. To 
pass from Eq. (23) to Eq. (24) the vanishing perturba- 
tions at Oa and Lb were enough to get rid of the bound- 
ary terms. The perturbations of velocity and acceleration 
are arbitrary at Oa and L b because there is no prescribed 
orbit either before Oa or after Lb, while the velocity and 
acceleration perturbations at L~and 0+ must vanish for 
a match with the histories. The condition Gi = de- 
fined by Eq. (26) yields the electromagnetic equations of 
motion with the Lienard-Wierchert-Lorentz force [3], as 
evaluated in the Appendix. The gradient for variations 
of trajectory 2 is obtained analogously, by discarding the 
integration over the past history of particle 2 and defin- 
ing a sub- functional 6*2(62, 62) obtained from the above Si 
(14) by exchanging particle indices. The Banach space 
for arbitrary variations of both trajectories respecting 
the EHBCs is the direct product N'^^\xi) (g) N^^\x2) = 
M'^\xi,X2) with the norm given by |6i, 62|iv(xi,x2) = 
sup II61II4 + sup [I62II4 + sup II61II4 + sup II62II4, which is 
the natural physical space of orbits satisfying the EHBCs. 

Action (13) is not defined for the superluminal ele- 
ments of N^^'>{xi) (which have a large norm |6i|jv(xi)) 
because it involves taking the square-root of a negative 
number. The above defined norm guarantees that suf- 
ficiently small neighborhoods of subluminal orbits con- 
tain only subluminal orbits (by Proposition 1), but the 
set of subluminal orbits is not closed because Cauchy 
sequences of subluminal orbits can converge to luminal 
orbits. Moreover, luminal orbits can have small neigh- 
borhoods containing superluminal orbits, for which again 
action (13) is not even defined. In the following we re- 
lax the parametrization-invariance and construct a sec- 
ond Poincare-invariant functional defined everywhere in 
A'^'^"'^) (a;i) and yielding the same electromagnetic equa- 
tions of motion. For superluminal trajectories the light- 
cone condition (4) can have an arbitrary number of zeros, 
and for these the double integration on the right-hand- 
side of Eq. (2) is generalized by extending formula (6) to 
all zeros of (4) in the integration interval (0, Xip), which 
prescribes a vanishing integral in the case of no solution 
in the interval. Our second functional is obtained by fur- 
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ther generalizing the kinetic terms, i.e., 

ft = - /"^(xi.Xi)frfAi-^/ ^(X2-X2rrfA2 



+ I M^j) • x.dXj. 





(27) 



The last term of action (27) is the double integral of ac- 
tion (13) written in a convenient form and extended to ar- 
bitrary orbits by evaluating A{xj) with Eq.(8) extended 
to all the zeros of the lightcone condition inside [0, Afe^]- 
Notice that the last term of (27) is still paramctrization- 
independent, unlike the generalized kinetic terms of (27) 
that are parametrization-invariant only ii p = 1/2. The 
Eulcr-Lagrange condition of extremum (25) applied to 
action (27) yields 

m~[{xi ■ Xir-'x>^] = E xUdkiA^ - d^iA% (28) 

"•M k=i 

where the partial derivative respect to the covariant com- 
ponents is defined by dki = d/dx^ and we expressed 
the Euler Lagrange condition (25) leaving the kinetic 
terms on the right-hand side. Equation (28) involves an 
anti-symmetric tensor on the left-hand side, so that the 
Minkowski scalar product of (28) with the four-velocity 
Xi yields zero on the left-hand side, i.e., 

i2p -1) d 



2p dX 



■pi 



= 0. 



(29) 



The Fokker-like action (13) has = 1/2 so that Eq. (29) 
holds trivially, but for p ^ 1/2 condition (29) implies 
that [xi - Xi) must be constant along the extremum orbit. 
Moreover, for p ^ 1/2 action (27) is no longer parame- 
ter independent, and Eq. (29) shows that the extremum 
condition of (27) is expressed in a parameter that along 
subluminal orbits is proportional to the proper-time pa- 
rameter (the constant of proportionality renormalizes the 
scalar mass of each particle). Property (29) divides the 
orbits in three invariant classes, as follows (a) if condi- 
tion (1) is positive at any orbital point, then it must 
be positive at all points of an extremal orbit, and (b) 
if the particle ever travels at the speed of light, then 
{xi ■ Xi) = everywhere along the extremal orbit, so 
that the particle travels at the speed of light everywhere 
and last (c) a superluminal orbit is superluminal every- 
where. By combining the three types of trajectory for 
each particle we can produce six different classes of or- 
bits, luminal 1-luminal 2, superluminal 1-luminal 2, and 
etc. . . We henceforth take p = 1 so that the kinetic in- 
tegrand (27) is analytic and more important action (27) 
is defined everywhere and Frechet-differentiable every- 
where in the Banach space A''(^^(a;i). The advantages of 
an action defined for trajectories violating (1) is that the 
ambient space of the functional (27) can be a complete 
normed linear space, even if we later decide that only 
subluminal orbits are interesting for physics. 



THE SECOND VARIATION 

Here we calculate the second variation about the 
C°°low-velocity-circular-orbit-cxtrcma (aj^.ajj) of large 
enough radii [9, 10] using either action (13) or (27). To 
calculate the first variation, in Ref. [11] we have ex- 
panded the delayed arguments of action (13), a method 
that becomes cumbersome for the second variation. Here 
we use a method motivated in the derivation of the low- 
vclocity-limit of the Fokkcr action [13], only that ours is 
not restricted to low-velocities and includes delay. Our 
method is equivalent to expanding the delayed arguments 
of action (13) but we use a shortcut equivalent to taking 
derivatives of the Dirac delta-function, as done in Ref. 
[13]. We start with a definition for the derivative of the 
right-hand-side of (6) and (7). The following proposi- 
tion justifies the formal manipulation of the S symbol 
inside integration-by-parts formulas as long as the inte- 
grand vanishes at the endpoints of the integration inter- 
val. To motivate our next definition we start from for- 
mulas (6) and (7) with trajectories given by a perturbed 
circular orbit, i.e.. d{\i, X2, e) = (ic(Ai, A2) +e?i(Ai, A2, e) 
with u(Ai,A2,£) given by a polynomial function of the 
bi(Ai) vanishing at the endpoints of [L21, L2f] according 
to (18). Definition 1 yields 

J ,5(d(Ai,A,£))/(Ai,A)dA (30) 



^ /(Ai,X 



^ (-^1 ) '^2''^ ) and the summation of Eq. 



where 2 Ja, = -g^ 

(30) is extended to all zeros (Ai, A2 ) of d{Xi, \2,e) with 
~ (7) 

A2 G [L2i,L2f] for any fixed Ai € [Lij,Lip]. The con- 
dition m(Ai, A2, e) = at the endpoints ensures that the 
lightcone condition is not perturbed at the endpoints, so 
that no zero A2 of d(Ai, A, e) leaves or enters the inter- 
val [L21, L2f] for small e. The implicit function theorem 

-(j) 

for d{Xi,X,s) = defines A2 as a function of e with 
derivative 



dX 



dd 
de 



(Ai,A 



0-) 
2 1 ' 



de 



g^(^Ai,A2 ,e) 



(31) 



The derivative of the right-hand-side of Eq. (30) respect 
to £ can be expressed with the help of (31) in the form 



9aJ./(Ai,A2)4/2JaJ | 
k ||2Ja,(Ai,A2)|| 



(32) 



where 2J\^ = and ds = |f(Ai,A2,e). Equa- 

tion (32) is formula (6) with /(Ai,A2) replaced by 
d\{f{Xi,X2)de/2J\^), an equality that justifies the use 
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of a formal derivative of the delta-function symbol as fol- 
lows 



de 



5(d(Ai,A,e))/(Ai,A)dA 



(33) 



L2I 



= J dJ{d{Xi,X,e)f{X,,X)dX 

= J dx{f{Xi,X2)d,/2Jx)6{d{Xi,X,e))dX 

where again 2Jx^ = — ^^(Ai,A2,e) and c?g = 
ff (Ai , A2 , e) vanishes at the integration limits. We hence- 
forth use (33) to define the formal derivatives of the delta 
symbol, stressing that there is no distributional limit in- 
volved but rather the above-defined operation. More- 
over, actions (27) and (13) depend on a double integral, 
i.e., either one of formulas (10). The derivative of the 
interaction / defined by Eq. (30) with an e-dependent 
d{Xi, X2,e) is given by either one of formulas 



de 



(34) 



Y^x dxA{xi-X2)de/2Jx,] ,-^(k) 

J dX2}^ |Ai 

L21 k \\2J\i\\ 

y dXi E ^^2Mi:±^}M^^\-xf 



I|2Ja. 



-^(Ai,A2,e), as long as d^ 



with 2J\. = — 7TT- 

ff(Ai,A2,£) = at the integration limits. Otherwise 
we might have to chose the line of Eq. (34) for which 
ff (Ai, A2,£) vanishes at the integration limits. 

For the second variation we vary both trajectories si- 
multaneously according to 

xi = x1+sbi, X2 = X2+eb2, (35) 
xi = xl+eb-i, X2 = X2+sb2, 

e G [0,1]. We require vanishing perturbations at the 
endpoints, 



61 (Oa) = 6i(L-) =0, 
62(0+) = 62(^5) =0, 



(36) 



and vanishing velocity and acceleration perturbations on 

the history side of each trajectory, i.e., at point O"*" of 
trajectory 2 and at point L~ of trajectory 1, 



bi(L-) = b2{L-) = 0, 
62(0+) = 62(0+) = 0, 



(37) 



so that the trajectories can be continued to a trajec- 
tory 61 = on the boundary segment (L~, L+) and 62 = 



on the boundary segment [O^ .O^). The quadratic 
integrand of the Taylor expansion involves products of 
variations at points connected by the lightcone condition 
(rather than variations at the same time as in the Kepler 
problem). We expand the action in a Taylor series up 
to the second order in e by using a directional derivative 
along the trajectory variation (35) of N'^'^^ {x'ljX'j). 
Once the circular orbit is an extremum, the first vari- 
ation vanishes so that Taylor's theorem gives the func- 
tional at £ = 1 as a sum of its value at e = plus the 
second- variation evaluated at some e € [0, 1]. 

The second variation of the first term on the right- 
hand-side of Eq. (13), representing the kinetic energy 



is 



1^ ■ 

/ 



. {x\-bxf-xfb\ 
2(i;J)3/2 



dAi 



(38) 



Formula (38_) is positive-definite, which is seen as 
follows; — If biis time-like, the positivity is given by 
the reverse-Schwartz inequality of time-like vectors men- 
tioned in the introduction, while for a space-like h\ Eq. 
(38) is a sum of positive terms. Since the interaction inte- 
gral is naturally expressed as a double integral times the 
Dirac delta-function we henceforth normalize all integrals 
to that form. To normalize Eq. (38) we simply add a 
dummy integration over dX2 multiplied by the integrating 
factor 2|| J;^2||c^A2 and use that 2||JA,||(5^rfA2 = 1, 
yielding 



A(2)ii'i = mi£2 y y dXidX2 

O- Oa 



(it-b,) 



xfb' 



(a;5)3/2 



(39) 

The symbol 5^ is an abbreviation for 6{d{Xi, X2,s)) as 
of definitions (6) and (7) while upper index c denotes 
the circular-orbit functions. Notice that the low-velocity- 
limit of IIJA2II ill particle-time parametrization is the 
spatial separation in light-cone, ri2, as defined by Eq. 
(12). To abbreviate notation we henceforth indicate the 
double-integral over both circular orbits of any integrand 
5(Ai,A2,e) times S{d{Xi, X2,s)) by J^g . For example 
the kinetic term Eq. (39) is abbreviated to 



/ 



{x'i-b.f 



xfbl 



(i;J)3/2 



(40) 



Next we calculate the second- variation of the interaction 
term by substituting variation (35) into the integrand 
If = S{\xi — X2^)x\ ■ X2 and expand in a Taylor series 
in £ using the above define rules for the formal derivative. 
The separation d(Ai,A2,£) is perturbed along variation 
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(35) to 



yields 



d{Xi,X2,e) = \xi - x-i^ 

+£'|6i-&2|', (41) 

so that the formal expansion of J/j becomes 

8n = 6d{\x\-x''^\^) 

- xl) ■ (bi - 62) + |bi - b2p]4 
+2[x52-(6i-62)]'<5L + 0(3), (42) 

where x^^ = {x\ — X2)- The bilinear product Xi ■ X2 is 
perturbed to 

Xi-X2=xl-X2 + 2s{xl • 61 + ±2 • ''2) + &i • b2 (43) 

Henceforth one or two primes over Sd denote respectively 
one or two formal derivatives as defined by formulas (33) 
and (34). The quadratic term of the Taylor expansion of 
If = S{\xi — cc2p)^i • X2 is obtained multiplying (42) by 
(43) and collecting the second order terms, yielding 



D 



A(2)7j, = e%-b2SD + e'^xl-x'2\bi-b2\'^5 
+2s^[xl2 ■ {bi - b2)]A{xi ■ X2)4 
+26^x1 ■ xl[xl^ ■ (6i - b2)]H'j,. (44) 

where A(a;i •a;2) = (cci -62 +^2 '^i) ^"^^ ^12 = (^1 "^2)- 
The first term on the right-hand side of Eq. (44) is al- 
ready in the normalized form of Eq. (39). Wc henceforth 
drop the factor of the second-order expansion. The 
second term on the right-hand side of Eq. (44) must be 
split in three monomials, {b\-2bi-b2 + bl)5 and 

the formal integration by parts to get rid of the Sj^ must 
treat each monomial differently cause formula (33) needs 
a vanishing perturbation at the endpoints;- For example 
the monomial xl ■ X2b\5]j must be dealt with according 
to the first line of (34), i.e.. 



J aAi^ 



d .{x'{-xl){bi-b2-b\)^ 



2||JaJ 



(46) 



J d\2^ 



d .(i;^i;^)(bi-52-b2) 



2|| Ja, 



Next integrating by parts on the second line of the right- 
hand-side of Eq. (44) yields 



J dXi^ 
c 

J dX2^ 



d ^A(d:|^x|)(x|2ibi), 

\\JxA\ ^ 

d ^/^{ xl ■ xl){xi2 ■ b2) . 



(47) 



where again A(£Ci ■ X2) = {x'{ •62 + 3^2 • bi) and x\2 = 
{xl — ajj). Last, the third line of the right-hand-side of 
Eq. (44) is transformed after two integrations by parts 
into 



1 /• 

2 J dXi^px, 
c 

4/ 



d xl ■ xl{xl2 ■ bi)^ 



]] 



9X1' \\Jx,\\ 
1 d xl- x%{x\2 ■ b2f 



(48) 



2i dX2'\\Jx^\\dX2 

c 

d . 1 



]] 



IkA.II 

d xl ■ xl{xl2 ■ &2)(a;i2 • 61) 



dXi'\\JxA\dX2' 



Henceforth we specify the 
particle-time parametrization, i.e. 



circular 
bi = 



]]• 



orbit 
(0,bi) 



adopting 
and bi = 



(0, bi) and bi ■ bj = hi ■ hj, where a dot between the 
vector parts henceforth denotes Cartesian product. The 
velocities along a limiting circular orbit of large radius 
are given by xl = {'^,v^) with 



-k: "^2 X 

^1 = Tr^^(*i)' 



(49) 



V2 = - 



mi 



V{t2). 



I 

Oa 



x\ ■ X2b\5]jdXi = 



(45) 



L+ 



Oa 



^\ X 9 ,Xi ■ ^2^1 ^ 



since the integrand on the left-hand side of Eq. (45) 

vanishes with bJ at L+and Oa (the EHBCs). Notice 
that the monomial with b2 does not vanish at L+and 
Oa- In that case the integration of choice would be over 
dX2. Using the above term- wise integration, the second 
term of the first line on the right-hand side of Eq. (44) 



In Eq. (49) ri2 is the constant separation in lightcone 
along the circular orbit, v(t) is the unit vector along the 
trajectory of particle 1 and M = mi + 1112 (the Kepler 
orbit is discussed in Ref. [11] ). The period of the circular 
orbit is given by Kepler's law 



r = 27r 



M 



mim2 



3/2 

''12 5 



(50) 



so that the lightcone separation ti = t2 i ri2 is a negli- 
gible fraction of the period for large ri2, i.e., the times 
in lightcone are almost equal {t\ — 12), the spatial posi- 
tions are almost in diametral opposition and the veloc- 
ities have nearly opposite directions. Using the above 
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circular orbit we calculate \\Jti\\ = \\Jt2\\ — ^12 and 
'^tilktill = '^t2ll'^f2ll = • ~ 1 in the limit of a 
large ri2. 

Theorem : The second variation about circular orbits 
of large enough radius is a strongly-positive quadratic 
form for trajectory variations satisfying (36) and (37). 

Proof: There arc three basic types of integrals of 
quadratic monomials in Eqs. (46), (47) and (48), 
namely (a) velocity-velocity, (b) position-position and 
(c) position-velocity. Notice that integrals of type 
/ (5_d(A • bi)(B • hj) can be re-expressed as an integral 
of a quadratic form of position and velocity variations 
only using (34). In the following we inspect each type 
of integral, finding that (a) is strongly-positive while (b) 
and (c) are dominated by (a) at large enough radii, as 
follows; 

(a) The velocity-velocity terms of the second-variation 
are 



A'-^^V = ^(miriab? -Fmzraiba + bi • bz), 



(51) 



which is strongly positive-definite at large separations, 
rriir 12 >> 1- (b) The dominant quadratic terms in the 
displacements are 



|bi - b2|^ ^ j- 3(n-bi-n.b2)^ ^g^) 



c c 
2(n-bi)(n-b2) 



+ 



+ 



/ 

c 

J 



12 



{V2 ■ hi){vi ■ b2) + {Vi ■ hi){V2 ■ b2) 



' 12 



where fi = (bi — b2)/ri2. Quadratic form (52) is not 
positive-definite, and in fact for bi = b2 = with 
n • t)= we have A^'^^R = -2mim2\\SR\\'^ /{Mrf2)- The 
first two lines on the right-hand side of Eq. (52) have a 
non-negative sum, while we can show using (49) that the 
last line is bounded, i.e.. 



miTO2 
" Mr?2 



ibiir+iib2in. 



(53) 



Lemma i;-For variations hi vanishing a± U = and 
U = Tff, it follows from the Fourier series hi = 
'^aksm{wkt/T^) that 



> 



(54) 



where bi = dhi/dU. In Eq. (54), is the time for the 
circular rotation to travel the angle 4> from Oa to L_, 
i.e., = ((/)/27r)T where T is the period as defined 
by Eq. (50). The equal sign in (54) holds iff the first 



Fourier mode alone is present, i.e., Ofe = for fc 7^ 1. 
The following inequality is true for arbitrary arcs of circle 
(f) <2-K but for simplicity we write it for EHBCs going a 
complete turn, Ta, = T, i.e.. 



ri2TOi||bi|P > 



mim2mi 
4Mrl2 



lb,; I 



(55) 



Using Eqs. (55) and (53) we can show that an arbi- 
trary fraction < / < 1 of the kinetic term (51) dom- 
inates the quadratic form (52) for sufficiently large ri2, 
i.e., /A(2)y > A(2)i?. 

(c) The quadratic terms involving position-velocity 
perturbations are also dominated by the kinetic terms, as 
follows; — Notice that the position- velocity terms coming 
from (46) integrate to zero, i.e.. 



1 



2ri2 



2ri2 



(bi • b2 - 2bi • bi) 



(b2 • bi - 2b2 • bs). 



(56) 



where we used the large-radius limits ||Jti|| = \\Jt2\\ — 
ri2 and i;^ • ~ 1 and moved ri2 outside of the inte- 
gration sign because it is constant along circular orbits. 
After integration over one parameter Eq. (56) reduces 
to the integration of an exact differential vanishing at 
the boundaries, so that (56) vanishes. The largest non- 
vanishing position-velocity terms come from (47) and 
(48), i.e.. 



A(2)vi? = (57) 

+ [ — A(vi ■ V2)(n • bi — n • b2) 
J ri2 



c 

■I 

c 



— A(i7i • V2){vi ■hi-V2- b2), 
ri2 



where n = (bi — b2)/ri2 and A{vi • V2) = (z^ • b2 -|- ^2 " 
bi). To show that the kinetic form (51) dominates the 

velocity-position quadratic terms for large enough ri2 we 
use inequality (55) to derive Lemma 2: 

j ri2{mi\\hi\f + mj\\hj\f) > J {ri2mi\\hi\f (58) 
c c 



+ 



J mim2mj||bj|p. 



> 



mim2mimj 



\hi\\\\h, 



where the last inequality is simply the completion of 
a binomial square. It can be verified with Eq. (49) 
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that the cocfScicnts of the monomials in the integrals 
of (57) are dominated by ^/r\!^, so that Lemma 2 as 
of (58) is enough for the kinetic terms to dominate all 
type (b) terms. To show that the second-variation is 
positive-definite we divide the kinetic energy (51) in 
three equal parts ;-The first third dominates the position- 
squared terms (52) for large enough ri2, as explained 
below Eq. (55), while the second third dominates the 
velocity-position terms by inequality (58). The last third 
is a non-degenerate positive-definite quadratic form of 
the velocities, so that the second variation about circular 
orbits of large enough radii is positive-definite, proving 
that circular orbits arc local minima. Moreover, the last 
third-part of Eq. (51) has all positive eigenvalues, so that 
the second variation is strongly positive. 

CONCLUSION AND DISCUSSIONS 

An important question is the existence of an extremiz- 
ing orbit for the functional (13) with arbitrary past data 
for particle 2 plus arbitrary future data for particle 1, 
i.e., the existence result for solutions of the mixed-type 
neutral-delay electromagnetic equations of motion with 
general boundaries. There are no existence results for 
the electromagnetic two-body problem apart from a few 
obtained for a one- dimensional motion with repulsive in- 
teraction [17], a qualitatively different and simpler case 
where the equations are not neutral but rather delay- 
only. For sufficiently small deformations of the circu- 
lar EHBCs preserving the boimdary lightcones Oa — O"*", 
Oa ~0~ and L~ —Lb and L+ —Lb, the second variation 
can be proved positive-definite with analogous methods. 
Moreover, on a subset O C A^'-^^ (a^J, ajj) of orbits satis- 
fying M > sup(||bi||4) -|- sup(||62||4) for some M we can 
reconstruct the perturbation using (18) and the one- 
sided conditions (37), a formula analogous to Eq. (19). 
For example for bi(Ai) we have 

6i(Ai) = /dAb / 6i(A„)rfA„, (59) 

Ac AiF 

from which it follows that sup(||6^'^'' ||) < A^^^*^ sup ||bi||4, 
with an analogous condition holding for 62 from the other 
side. Conditions (37) can be used to show that the 
perturbations inside O are equicontinuous and uniformly 
bounded, so that by the Arzela-Ascoli theorem the set 
6 is compact. If the second variation is positive-definite, 
the functional is bounded from below on the compact set 
and assumes its minimum inside O. We conjecture 
that this point of minimum is an interior point of the 
compact set. That granted, the minimum has a whole 
neighborhood inside 6, so that Eq. (24) holds for arbi- 
trary bfc and the gradients must vanish at the minimum. 
Conditions Gk = with Gk defined by Eq. (26) are 
the state-dependent neutral-delay equations of motion. 



so that this would be an existence result for the state- 
dependent neutral-delay equations. This result would be 
the analogous of the "Kurtzweil small delays don't mat- 
ter theorem" for global trajectories of DDE's on com- 
pact sets[16]. The uniqueness theory also differs from 
the case of Ref. [21], and here one should again start 
from the case of slightly perturbed circular boundaries, a 
case where the equations of motion are approximated by 
neutral-delay-equations with constant advance and delay. 
For circular orbits of intermediate radius some inspection 
suggests the minimum should become a saddle in a bi- 
furcation at a finite 0(1) radius in our unit system, i.e., 
of the order of the classical electronic radius. 

The existence proof is much harder for the solenoidal 
orbits discussed in the appendix because of the denom- 
inators. For solenoidal orbits with a fast velocity the 
functional might have a maximum as suggested by the 
kinetic term. 

A useful generalization of our second functional is for 
orbits defined on a Sobolev space Hq with derivatives de- 
fined almost everywhere. For that we need to generalize 
the lightcone condition to arbitrary trajectories and to 
generalize Eq. (6) to a sum over all zeros of the lightcone 
condition. The fact that Eq. (6) is further integrated 
over the other orbital parameter to make Eq. (10) com- 
pensates for the extra zeros gained by changing the tra- 
jectories on a set of zero measure, so that the functional 
can be defined on Hq. This generalization could be useful 
in proving existence for the case of general boundaries. 

Another question of interest regards the search for pe- 
riodic orbits and the possibility to restrict the variational 
method to the family of periodic orbits satisfying the 
EHBCs. The reduction is possible to a sub-family of pe- 
riodic orbits by identifying the spatial components of Oa 
with those of for trajectory 1 and the spatial com- 
ponents of O- with those of L^ for trajectory 2, which 
must be the case along a periodic orbit. The orbital vari- 
ation inside the family of periodic orbits must preserve 
the history segment of each trajectory, as illustrated in 
red in Fig. 3.1, which is a sub- family of the family of 
periodic orbits. Last, it is possible to extremize the func- 
tionals directly in the space of G^ orbits without even 
respecting the former sub-family conditions. The condi- 
tions for an extremum with these most general variations 
are no longer the electromagnetic equations of motion but 
rather the overdetcrmincd equations obtained by vanish- 
ing both linear terms on the right-hand side of Eq. (23) 
separately. 
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APPENDIX: PHYSICS OF THE FOKKER 
ACTION 

Here we evaluate the gradient (26) explicitly and dis- 
cuss the physics of the two-body problem using proper- 
time parametrization for the trajectories. The velocity 
respect to proper-time can be expressed either in the form 
(21) with hi = 7j, i.e., 



or in the form 



where 



dTi 



7i = > 0. 

dTi 



(60) 



(61) 



(62) 



According to Eq.(l) the velocity respect to proper-time 
along physical orbits, Vi = dxi/dri, satisfies 



{Vi ■ Vi) = 1, 



(63) 



which can be solved for 7^ using either Eq.(61) or Eq.(60), 
yielding 



7i 




(l-|Kin-V2. (64) 



positive number ti - t2- = 7'12- ^ 11^2(^2-) - ri(ti)\\. 
The same Eq.(5) with the plus sign (for the advanced 
cone) yields the negative number ti — t2+ = —ri2+ = 
— 11^2(^2+) — n.(ti)||- The sign of the time-component is 
explicitly indicated by the plus or minus on the first en- 
try of Eq.(65). The half-Jacobian (11) with proper-time 
parametrization is calculated using Eqs. (60) with index 
1 replaced by 2 and definition (65), i.e.. 



= (a;i2 • X2±) = T72»'i2±(l ± n ■ V2)±, 



(67) 



where overdot represents derivative respect to proper- 
time of particle 2. Notice that can become singular 
when the particle moves near the speed of light. 

The partial derivative of with respect to xi in 
Eq.(26) along the fixed trajectory of particle 2 is given 
by formula (16) with dtf replaced by drf. Since we are 
operating with proper-time wc can set \/xi ■ Xi = 1 in 
the first denominator on the right- hand- side of Eq. (26), 
yielding 



miXi 

d ^ ±2- 



(68) 



d 



axi 



T2 

Xi ■ X2 



X2+ 

2JX 



Xi ■ X2+ ~ 



2Jr 



2J. 



+ 



where we took out the modulus sign using that is 
negative and is positive. The derivative respect to 
Ti on the right-hand side of Eq.(68) also acts on the 
arguments r J since these are functions of ri by the light- 
cone conditions 



\xi{ti) - X2{t 



0. 



(69) 



To evaluate the derivative of the retarded/advanced 
proper-time with respect to t\ we take the differential 
of the light-cone condition (69), i.e.. 



Notice that ||^|| is unbounded and becomes arbitrar- 
ily large when the time- velocity approaches the speed of 
light. Using condition (3) to solve for the Euclidean norm 
of the spatial separation, we can express the separation 
vector X12 = {xi — X2) as 



X12 = {Tn2,n2n'^ 



(65) 



where ri2± is the distance in light-cone and n='= is defined 
by 



.± ri(ii) - f2(i2±) 

n = — ^ — ^ 

ri2± 



(66) 



a unitary Euclidean three-vector. Notice that the spatial 
distance in light-cone ri2± = ||?^2(i2±)— ?l(ii)|| is adiffer- 
ent function for each light-cone. The time-component of 
X12 = {xi — X2) is simply ti — t2 and evaluated with the 
negative sign of Eq.(5) (for the retarded cone) yields the 



2(xi2 • xi)±dTi - 2{xi2 ■ X2)±dTf = 0, (70) 



which yields 



drf {xi2 ■ Xi)- 



dri {xi2 ■ X2)- 



(71) 



where Xx2± = {xx — X2±)- Formula (71) is vaUd for both 
the retarded and the advanced lightcones. The same sep- 
aration X12 appears on both numerator and denominator 
on the right-hand-side of Eq.(71), so that the plus or mi- 
nus sign of Eq.(65) cancels out and the derivative (71) is 
always positive as it should be. Wc can use the two signs 
of Eq.(71) to calculate the derivative of the most retarded 
argument with respect to the most advanced argument 
by the chain rule 



dT2- _ , dT2- ., dn . 
dT2+ ~ ^ dn '^dT2+'' 



(72) 
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where = ( ) ^ and Eq.(72) is a non- negative 

rate because it is a product of two positive factors. 
The fact that the retarded and the advanced arguments 
have non-negative rates ensures the continuabiHty of any 
piecewise-continuous solution at a breaking point [22]. 
Therefore the usual mechanism for a neutral-delay equa- 
tion to loose its piece-wise continuous solution at a break- 
ing point is absent and the neutral equations of electro- 
dynamics never loose solutions for this reason. 

Using Eq.(71), the second line on the right-hand side 
of Eq.(68) evaluates to 

' -X2 - (xi2 • a;i)a2j- (73) 



2{X12 ■ X2)- dTi 



1 .d{Xi2 ■±2) . , . . -, 

' -X2 - {Xi2 ■ Xi)a2\ + . 



2{xi2 ■ X2)\ dn 

where the lower index ± after the bracket indicates eval- 
uation in the advanced/retarded light-cone respectively 
and a2± = d^X2/dT2± denotes the acceleration of par- 
ticle 2 respect to proper time in the advanced/retarded 
light-cone respectively. Last, on the third line of the 
right-hand side of Eq. (68) the partial derivative respect 
to xi acts on xi and also on quantities of particle 2 by 
the rule 



d 
dxi 



dT2 d 

dXi dT2 



(74) 



with given by Eq.(16). The manipulations are sim- 
ple and the third line on the right-hand side of Eq.(68) 
becomes 



{Xi-a2)-Xi2- {Xi-X2)-X2- 



+ 



2(a;i2-i;2)^ 2(a;i2 • ±2)^ 
{xi ■ X2)-[l - (a2 • a;i2)-]a;i2- 

2(a;i2 • ±2)- 
(xi ■ a2)+a;i2+ _^ (cfci • X2)+X2+ 
2{xi2-X2)\ 2(a;i2-i;2)+ 
(Ai ■ a;2)+[l - (02 ■ a;i2)-F]a;i2+ 

2(£Ci2 • X2)\ 



(75) 



(76) 



Using Eqs.(68), (73) and (75) we can express the gradient 
as 



Gi = mixi - ^F} - ^F2 , 



(77) 



where 



^± = {Xl2-V2)± 



+ 



± [{xi2 ■ vi)a2 - {vi ■ a2)xi2]± (78) 

[(a;i2 • ^1)1^2 - {vi ■ t^2)a;i2]±, 



2P2 

(1 - X12 ■ a2)± 



where/9^ = || and V2± = dx2/dT2±- 

The condition Gi = yields a familiar Newtonian-like 
equation of motion with the Lorentz-force of the other 



particle as a semi-sum of advanced/retarded Lienard- 

Wicchcrt fields, i.e., 

"^1^ = ^F2{xi,Vi,X2iT2),V2iT2),ai2{T2)) (79) 

+ ^^2"^(a'l> ^1: a;2(T^), W2(t+), a2(T+)), 
TO2^ = ^F^{x2,V2,Xi{T^),Vi{T^),ai{T:[)) (80) 

aT2 ^ 

+ ^F^{X2,V2,Xi{Tt),Vi{Tt),ai{Tt)), 

where Xi,Vi, Oj are respectively the position, velocity and 
acceleration of particle i with respect to proper-time Tj. 
In Eqs. (79) and (80) the forces depend respec- 
tively on the other particle's retarded/ advanced posi- 
tion, velocity and acceleration, as well as on the object- 
particle's present position and velocity. Moreover the 
retarded/advanced points are implicitly defined by Eq. 
(5), so that Eqs. (79) and (80) are neutral-delay equa- 
tions of mixed-type with implicit state-dependent delay. 
The forces F^ of Eqs. (79) and (80) are the Lorentz 
force of the Lienard-Wierchert fields of standard electro- 
dynamics textbooks [6, 23]. Notice that each line of Eq. 
(78) is orthogonal to Vi, so that it follows from Eq. (79) 
that 



{vi ■ vi) = 0, 



(81) 



in agreement with Eq. (29) for p 7^ 1/2 and for p = 1/2 
Eq.(81) is the definition of proper-time parametrization. 
Condition (63) can be solved for the time-velocity with 
Eq. (64), thereby reducing the dynamics to the spatial 
components of Eqs. (79) and (80). Last, in the follow- 
ing we discuss the denominators of Eqs. (79) and (80). 
The Lienard-Wierchert force (78) involves denominators 
of type (67), which become singular when the other par- 
ticle travels near the speed of light. This is illustrated 
expressing the vector-part of the equation of motion (79) 
using particle-time parametrization and using only the 
leading singular term of force (78), i.e., 



-1 



zTTo-n X (n X dj) + ... (82) 



d , niiVi _ ^ 

Notice that the left-hand-side of Eq. (82) becomes sin- 
gular when particle i travels near the speed of light 
while the right-hand side of Eq. (82) becomes singu- 
lar when particle j travels near the speed of light in ei- 
ther the past /future lightcone points, so that if the two 
motions synchronize a solenoidal orbit with a fast ve- 
locity could exist, as suggested in Ref. [11]. Surpris- 
ingly this non-trivial motion does not require large total 
momenta, as follows; — Action (13) is invariant by the 
Lorentz group if one also moves the boundary-condition- 
segments with the group element (the red segments of 
Fig. 3.1). Noether's theorem [13] applies to action (13) 
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in a way completely analogous to the formal derivation of 
Scliilci[10], as explained in Ref. [13], yielding invariants 
defined by finite integrals, i.e.. 



P5'(n)+P^(r2) + 



(83) 



no numerical integrator available to integrate such non- 
trivial state-dependent neutral-delay equations, and we 
hope this work contributes to the construction of such 
integrator. 



L+ T2 



-2 



ri O- 
Ti Lb 



+2 



(CC^* - X^)5'{\xi - X2\^)X1 ■ X2dTidT2 

- a;2)i5'(|a;i - X2\'^)xi ■ X2dTidT2- 



Oa 1-2 



[3] 



and 



[4] 

(r?pf-r?p?)|.,+(r^pf-rfp^)|., [5] 

-2 y y {X^X^ - xfx^)5'{\xi - X2\'^)X1 ■ X2dTidT2 [7] 

- f [8] 

ri Lb 

+2 / [ {X^X^ - X^X^)6'{\Xi - X2f)Xi ■ X2dTidT2 [9] 

[10] 

[11] 

[12] 



Oa 1-2 



+ I I (A^if - xfx^)d'{\xi - X2\^)dTidT2 



Tl O- 
Ti Le 



(i^ccf - x"^x'^)6'{\xi - X2{^)dTidT2. (84) 



[13] 
[14] 



Oa r2 



where 



P'2 



Vl-I^il' 



'2+ 



^''12(1 ~ ^ ■ ^2-) 2rj'"2(l — n • U2+) 



(85) 



(86) 



'1+ 



2r2i(l — n • iTi_) 2r2i(l — n • tTi+) 

Notice that pi and p2 as defined by Eqs. (85) and (86) 
can be small even at fast velocities, so that a solenoidal 
motion with a stiff gyration near the speed of light is 
possible with finite and small mechanical momenta (83), 
as illustrated in Fig. 6.1. The solenoidal orbits of Ref. 
[11] were estimated to have a velocity near the speed 
of light, and the physical interest stems from the fact 
that these can be found in the physical region of small 4- 
momentum and angular-momentum. At present there is 
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